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Abstract. We present ISO photometric and spectro- 
scopic observations of a sample of 57 bright Asymptotic 
Giant Branch stars and red supergiants in the Large Mag- 
ellanic Cloud, selected on the basis of IRAS colours indica- 
tive of high mass-loss rates. PHOT-P and PHOT-C pho- 
tometry at 12, 25 and 60 /xm and CAM photometry at 
12 /nm are used in combination with quasi-simultaneous 
ground-based near-IR photometry to construct colour- 
colour diagrams for all stars in our sample. PHOT-S and 
CAM-CVF spectra in the 3 to 14 /im region are pre- 
sented for 23 stars. From the colour-colour diagrams and 
the spectra, we establish the chemical types of the dust 
around 49 stars in this sample. Many stars have carbon- 
rich dust. The most luminous carbon star in the Mag- 
ellanic Clouds has also a (minor) oxygen-rich component. 
OH/IR stars have silicate absorption with emission wings. 
The unique dataset presented here allows a detailed study 
of a representative sample of thermal-pulsing AGB stars 
with well-determined luminosities. 

Key words: Stars: carbon - circumstellar matter - Stars: 
mass loss - Stars: AGB and post-AGB - Magellanic 
Clouds - Infrared: stars 
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1. Introduction 

One of the least expected achievements of the Infra-Red 
Astronomical Satellite (IRAS; Neugebauer et al. 1984) was 
the detection of a large number of mid-IR point sources in 
the Large Magellanic Cloud (LMC) just above its limits of 
sensitivity (IRAS Point Source Catalogue; Schwering & Is- 
rael 1990). Many of these are candidates for intermediate- 
mass stars at the tip of the Asymptotic Giant Branch 
(AGB). Their lives drawing to a close, these stars are shed- 
ding their stellar mantles at rates of up to 10 -4 M Q yr _1 . 
Their dusty circumstellar envelopes (CSEs) obscure the 
optical light from the star and become very bright IR ob- 
jects. The details of the evolution and mass loss of AGB 
stars are poorly understood. The study of galactic samples 
of AGB stars is severely hampered by the difficulty to de- 
termine accurate distances to stars in the Milky Way. The 
distance to the LMC, however, is well known and hence 
luminosities and mass-loss rates of AGB stars in the LMC 
may be determined with a high degree of accuracy. 

Early explorations of the IRAS data in combination 
with ground-based near-IR observations resulted in the 
first identifications of mid-IR sources in the LMC with 
obscured AGB stars (Reid et al. 1990; Wood et al. 1992). 
We have successfully increased the sample of known AGB 
counterparts of IRAS sources in the LMC from a dozen 
to more than 50 stars (Loup et al. 1997; Zijlstra et al. 
1996; van Loon et al. 1997, 1998a: Papers I to IV). We 
attempted to classify their photospheres and CSEs as 
oxygen- or carbon-dominated, but for the majority of the 
stars this could not be done conclusively. There remained 
therefore considerable uncertainty about the luminosity 
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distribution of the obscured carbon stars. This informa- 
tion is important for testing current understanding of the 
evolution of AGB stars, including dredge-up of carbon and 
nuclear burning at the bottom of the convective mantle 
(Hot Bottom Burning, HBB). 

57 obscured AGB stars and a few red supergiants 
(RSGs) in the LMC were selected for Guaranteed Time 
and follow-up Open Time observations with the Infrared 
Space Observatory (ISO; Kessler et al. 1996). The goals 
were to obtain photometry at 12, 25 and 60 /im and 
to spectroscopically determine the chemical types of the 
CSEs. The photometry, which covers the entire spectral 
energy distributions (SEDs), can be modelled and used to 
derive accurate luminosities and mass-loss rates. In this 
paper we present the ISO data and classify sources as 
oxygen- or carbon-rich. 

2. Source selection 

The sources observed with ISO were selected from the 
lists presented in Paper I, where all IRAS candidate 
AGB stars in the MCs are listed. We selected 30 in- 
frared AGB stars or RSGs without optical counterparts 
from their Table 2. These objects should represent the 
brightest, most obscured AGB stars. Four objects from 
this table were excluded because of their red IRAS colours 
(SWSi2^2.5): LI-LMC528, 861, 1137 and 1341. We also 
selected 8 sources from the optically known M and C stars 
with IRAS counterparts in Table 1 of Paper I. These in- 
clude well known Harvard Variables as well as the opti- 
cally thick source IRAS04553-6825 (LI-LMC181, WOH 
G64; Elias et al. 1986; Wood et al. 1986). Two uniden- 
tified IRAS sources from Table 4 of Paper I have been 
included in the present sample. LI-LMC203 is near an 
Ml. 5 star (HV12501), but there is also an A3 lab su- 
pergiant (Sk— 69-39a) close to the IRAS position. For LI- 
LMC1795 we found a bright R-band counterpart (Paper 
II). Finally one source from Table 7 of Paper I was in- 
cluded (LI-LMC1130). Although listed in Paper I as a 
foreground star, it was included here in an attempt to es- 
tablish whether this is true. For these last three stars the 
higher spatial resolution ISO observations at 12 /im allow 
a better identification of the source with one of the pos- 
sible counterparts found near the IRAS position. The 41 
IRAS sources included in this study are listed in Table 
1, with the most common names for these objects, their 
coordinates (J2000) and some references. The coordinates 
for the pointings of the ISO observations were taken from 
the SIMBAD astronomical database in 1994. 

The selection of IRAS detected AGB stars gives a sam- 
ple that is severely biased towards very luminous stars (in- 
cluding supergiants). We therefore also included 16 non- 
IRAS stars. These were mostly taken from Wood et al. 
(1983, 1985), Reid et al. (1990) and Hughes (1989). Seven 
of these objects are classified as C stars from optical spec- 
tra or near-IR colours. Six objects are classified as M or 



S stars and for three objects no classification is avail- 
able. This group of non-IRAS sources also includes the 
RCB-like variable HV2379 (Bessell k Wood 1983). These 
sources are listed in Table 2. 

3. IRAS data 

Here the IRAS data are discussed, for later comparison 
with the ISO photometry. Data at 12, 25, 60 and 100 /im 
was retrieved from the IRAS data base server in Gronin- 
gen[] (Assendorp et al. 1995). The Groningen Gipsy data 
analysis software was used to measure the flux density 
from a trace through the position of the star (Gipsy com- 
mand SCAN AID). For the 60 and 100 fim data, 2x2 
square degree maps were created with 0.5' pixels to find 
point sources coincident with the positions of the stars. 
The 12 and 25 /im flux densities have a 1-a error of a few 
per cent, with a minimum error of ~ 0.01 Jy. The 60 and 
100 fim flux densities are much less certain, and it is also 
more difficult to assess reliable error estimates: 10% would 
be a typical error. The faintest 60 /im sources that IRAS 
detected were assigned F 60 = 0.1 Jy. Only one source was 
well detected at 100 /im. The flux densities are listed in 
Table 3. When it is not certain that the measured flux 
density is physically related to the star of interest it is 
marked with a colon. 

All of our sources that are in the IRAS-PSC, plus 
HV5870 (=LI-LMC1145) and TRM72 (=LI-LMC578) 
that are in Schwering & Israel (1990), were recovered 
with good flux density determinations at 12 /im. Reli- 
able 12 /im flux densities could also be determined for 
IRAS05128-6455 and 05289-6617, below their upper lim- 
its as listed in the PSC. Neither in the PSC, nor in Schw- 
ering & Israel (1990), are HV12070, HV2379, HV2446, 
TRM45, and TRM88 secure detections. Detection is not 
certain for WBP14 and the SHV sources for which flux 
density estimates are listed. The 12 /im flux densities of 
the GRV source and four SHV sources are upper limits. 
IRAS05506-7053 looks extended or multiple. 

At 25 /im detections seem a little more reliable than 
at 12 /im, at a given flux density. Rather surprisingly, the 
detection limit at 25 /im is at least as faint as at 12 /im; 
sources with F25 ~ 0.02 Jy could be found (see also Reid 
et al. 1990). This is, however, only possible because the 
positions of the stars are known. For SP77 30—6 and all 
eight (other) IRAS sources the PSC lists only upper limits 
of F 25 < 0.25 Jy. SHV0502469-692418 and WBP14 were 
the only sources that were (tentatively) detected at 12 /im 
but not at 25 /im. Their flux densities are probably below 

1 The IRAS data base server of the Space Research Organisa- 
tion of the Netherlands (SRON) and the Dutch Expertise Cen- 
tre for Astronomical Data Processing is funded by the Nether- 
lands Organisation for Scientific Research (NWO). The IRAS 
data base server project was also partly funded through the 
Air Force Office of Scientific Research, grants AFOSR 86-0140 
and AFOSR 89-0320. 
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Table 1. IRAS detected stars observed with ISO: names (LI stands for LI-LMC (Schwering & Israel 1990), TRM is 
from Rcid ct al. (1990), HV is from Payne-Gaposchkin (1971), SP is from Sanduleak & Philip (1977) and WOH is from 
Westerlund et al. (1981); sources will be referenced hereafter by their bold-faced names), ISO pointing coordinates 
(J2000), and references: 1: Hodge & Wright (1969); 2: Eggcn (1971); 3: Wright & Hodge (1971); 4: Dachs (1972); 5: 
Sandage & Tammann (1974); 6: Glass (1979); 7: Humphreys (1979); 8: Blanco et al. (1980); 9: Feast et al. (1980); 10: 
Bessell & Wood (1983); 11: Wood et al. (1983); 12: Rebcirot et al. (1983); 13: Prevot et al. (1985); 14: Elias et al. 
(1985); 15: Wood et al. (1985); 16: Elias et al. (1986); 17: Wood et al. (1986); 18: Rcid ct al. (1988); 19: Rcid (1989); 
20: Hughes (1989); 21: Hughes & Wood (1990); 22: Reid et al. (1990); 23: Hughes et al. (1991); 24: Wood et al. (1992); 
25: Roche et al. (1993); 26: Groenewegen et al. (1995); 27: Zijlstra et al. (1996); 28: Ritossa et al. (1996); 29: van Loon 
et al. (1996); 30: van Loon ct al. (1997); 31: Loup et al. (1997); 32: Oestreicher (1997); 33: van Loon ct al. (1998a); 
34: Groenewegen & Blommaert (1998); 35: van Loon et al. (1998b); 36: van Loon et al. (1999) 



LI 


IRAS 


TRM 


HV 


RA (2000) 


Decl (2000) 


Other names 


References 


IRAS detected stars 


1825 


04286 6937 


— 


— 


04 28 30.3 


-69 30 49 


— 


27,31,33 


1844 


04374 6831 


_ 


— 


04 37 22.8 


-68 25 03 


— 


27,31,33 


4 


04407 7000 


— 


— 


04 40 28.4 


-69 55 13 


— 


27,31,33 


57 


04496 6958 


— 


— 


04 49 18.6 


-69 53 14 


— 


27,31,33,34,36 


60 


04498-6842 


— 


— 


04 49 41.4 


-68 37 50 


— 


27,31,33,36 


77 


04509 6922 


— 


— 


04 50 40.2 


-69 17 33 


— 


24,27,28,33,36 


92 


04516 6902 


— 


— 


04 51 28.4 


-68 57 53 


- 


24,27,33 


121 


04530 6916 


— 


— 


04 52 45.3 


-69 11 53 


— 


24,27,28 


141 


04539 6821 


— 


— 


04 53 46.3 


-68 16 12 


— 


27,31,33 


153 


04544-6849 


— 


— 


04 54 14.4 


-68 44 13 


SP77 30-6, WOH SG66 


12,13,20,21,27,31 


159 


04545 7000 


— 


— 


04 54 09.8 


-69 56 00 


— 


24,27 


181 


04553-6825 


— 


— 


04 55 10.1 


-68 20 35 


WOH G64 


16,17,19,24,25,27,29,31,33,35,36 


198 


04557 6753 


— 


— 


04 55 38.9 


-67 49 10 


— 


27,31,33 


203 


04559-6931 


_ 


12501 


04 55 41.6 


-69 26 25 


SP77 31-20, WOH SG097 


11,12,13,20,22,27,32,33 


297 


05003 6712 


— 


— 


05 00 18.9 


-67 08 02 


— 


27,30,31,33,36 


310 


05009 6616 


— 


— 


05 01 03.8 


-66 12 40 


— 


27,31,33,36 


383 


05042-6720 


48 


888 


05 04 14.3 


-67 16 17 


SP77 29-33, WOH SG140 


5,6,7,11,14,18,22,31,32 


570 


05112 6755 


4 




05 11 10.1 


-67 52 17 




22,27,31,33,36 


571 


05113 6739 


24 




05 11 13.7 


-67 36 35 




22,27,31,33 


578 




72 




05 11 41.2 


-66 51 12 




22,27,31,33 


612 


05128-6728 


43 


2360 


05 12 46.4 


-67 19 37 


SP77 37-24, WOH SG193 


3,5,6,7,11,14,18,22,31 


1880 


05128-6455 






05 13 04.6 


-64 51 40 




27,31,33,36 


663 


05148-6730 


36 


916 


05 14 49.9 


-67 27 19 


SP77 37-35, WOH SG204 


1,2,4,6,7,11,18,22,31,32 


793 


05190 6748 


20 




05 18 56.7 


-67 45 06 




22,27,31,33 






88 




05 20 20.9 


-66 36 00 




22,27,31,33,34,36 






45 




05 28 16.3 


-67 20 55 




22,27,31,33 


1157 


05295 7121 






05 28 40.8 


-71 19 13 




27,31 


1130 


05289 6617 


99 




05 29 02.6 


-66 15 31 




22,27,31 


1145 






5870 


05 29 03.5 


-69 06 47 


SP77 47-17, WOH SG331 


9,11,20,31 


1153 


05294 7104 






05 28 47.8 


-71 02 29 




24,27,31 


1164 


05298 6957 






05 29 24.5 


-69 55 14 




24,27,31,36 


1177 


05300 6651 


79 




05 30 04.2 


-66 49 23 




22,27,31,36 


1238 


05316-6604 


101 




05 31 45.9 


-66 03 51 


WOH SG374 


22,27,31 


1281 


05327-6757 


5 


996 


05 32 36.0 


-67 55 08 


SP77 46-59, WOH SG388 


7,11,17,18,22,31 


1286 


05329 6708 


60 




05 32 52.5 


-67 06 25 




17,22,24,26,27,31,33,36 


1345 


05348-7024 






05 34 16.1 


-70 22 53 




27,31,33 


1382 


05360 6648 


77 




05 36 03.3 


-66 46 47 




22,27,31,33 


1506 


05402 6956 






05 39 44.6 


-69 55 18 




24,27 


1756 


05506 7053 






05 50 09.1 


-70 53 12 




27,31,33 


1790 


05558 7000 






05 55 20.8 


-70 00 05 




27,31 


1795 


05568 6753 






05 56 38.7 


-67 53 39 




27,31 
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Table 2. The list of program stars without IRAS counterpart. The references are as in Table 1. SHV is from Hughes 
(1989), BMB is from Blanco et al. (1980), WBP is from Wood et al. (1985) and GRV is from Reid et al. (1988). 



LI 


IRAS 


TRM 


HV 


RA (2000) 


Decl (2000) 


Other names 


References 


non-IRAS sources classified as C stars 










04 53 59 7 


-67 45 47 


SHV0454030 675031 


20,21 










05 02 28.7 


-69 20 10 


SHV0502469-692418 


20,21,23 








2379 


05 14 46.3 


-67 55 47 




3,10,11,20 










05 20 46.8 


-69 01 25 


SHV0521050-690415, BCB-R046 


8,20,21,23 










05 25 30.6 


-70 09 13 


SHV0526001-701142 


20,21 










05 26 17.4 


-69 08 07 


WBP14 


15 










05 35 11.4 


-70 22 46 


SHV0535442 702433 


20,21 










non-IRAS 


sources classified as M or S stars 










2446 


05 20 01.5 


-67 34 43 


WOH G274, GRV0520-6737 


11,18 










05 21 33.1 


-70 09 56 


SHV0522023 701242 


20,21 










05 21 40.5 


-70 22 31 


SHV0522118-702517 


20,21 










05 24 31.3 


-69 43 25 


SHV0524565 694559 


20,21 










05 30 00.3 


-70 20 06 


SHV0530323 702216 


20,21,23 








12070 


05 52 27.8 


-69 14 12 


WOH SG515 


9,11 


non-IRAS sources without spectral classification 










05 00 11.2 


-68 12 48 


SHV0500193 681706 


20,21 










05 00 13.5 


-68 24 56 


SHV0500233 682914 


20,21 










05 19 41.8 


-66 57 50 


GRV0519 6700 
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the limit of detection, F25 < 0.01 Jy, if their colours are 
rather blue. 

At 60 /im, IRAS05298-6957 is a bright, small but ex- 
tended source about 10' in diameter, with Fqo ~ 2 Jy. 
No point source could be distinguished on top of this 
emission, that is probably associated with the small clus- 
ter of which IRAS05298-6957 is a member (Paper IV). 
Flux densities are listed for two dozen sources, but it is 
not sure how many among these are real detections and 
how many are spurious. The only 60 /im detections in the 
PSC are IRAS04516-6902 (0.80 ± 0.19 Jy), 04530-6916 
(20.51 ±1.85 Jy) and 05112-6755 (0.91 ±0.11 Jy), all con- 
sistent with our estimates. More stringent upper limits are 
put on the 60 /im flux densities of the other sources. 

At 100 /im sources may be detected as faint as a few Jy. 
The only detection, however, is the brightest far-IR source 
in our sample, IRAS04530— 6916, which we measured at 
-Fioo = 28 Jy. This is consistent with the PSC upper limit 
of 46. 17 Jy. 

The new flux density estimates can be compared with 
the literature values from the PSC or Schwering & Is- 
rael (1990) (Fig. 1). On average, the new flux densities 
are only a few per cent fainter than the values from the 
literature. Flux densities F25S1O.2 Jy may have been over- 
estimated in the past. HV12501 with F lev /F\i t = 0.56 and 
0.32 at 12 and 25 /im, respectively, and IRAS05506-7053 
with F rcv /Fnt = 0.67 and 0.42 at 12 and 25 /im, respec- 
tively, are the most extreme examples of this. Schwer- 
ing & Israel (1990) over-estimated the 25 /im flux den- 
sity of TRM72 (F rcv /F lit = 0.55), but under-estimated 
the 12 /im flux density of HV5870 (F rcv /F lit = 2.00). 
The other two flux densities which are obviously under- 
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Fig. 1. A comparison of the estimates of IRAS flux den- 
sities given here and values from the IRAS Point Source 
Catalogue and Schwering & Israel (1990), at 12 (solid sym- 
bols) and 25 /im (open symbols). 



estimated are for SP77 30-6 at 12 /mi (F rev /Fi it = 1.53) 
and IRAS04286-6937 at 25 /im (F rcv /F m = 1.67). 

IRAS counterparts not listed in IRAS-based catalogues 
may still be found in the original IRAS data. This is be- 
cause manual extraction and measuring of the data is a 
more sophisticated technique than the automatic tech- 
niques that created the existing catalogues. In particu- 
lar, manual flux-density determination enables the back- 
ground flux levels to be estimated and subtracted bet- 
ter, yielding more reliable photometry. The only sources 
from our ISO sample that could not be detected in the 
IRAS data at either 12, 25 or 60 /im are GRV0519-6700, 
SHV0454030-675031 and SHV0524565-694559. 
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Table 3. Revised IRAS 12, 25, 60 and 100 /an photometry 
(in Jy), accompanied by a colon if questionable. 
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4. ISO observations 

The programme stars were observed with ISO at 12, 25 
and 60 /an (chopped measurements) and with PHOT-S 
as part of a Guaranteed Time programme under propos- 
als NTMCAGB1 and NTMCAGB2, and at 60 ^m with 
mapping mode and CAM-CVF as part of an Open Time 
programme under proposal LMCSPECT. 

The 12 fim photometry was mostly obtained using 
the ISOCAM instrument (Cesarsky et al. 1996) in staring 
mode with a 3" pixel field of view in beam LW-s and using 
the LW10 filter (~ IRAS 12 /an). We did 25 exposures of 
each 2.1s duration, after a number of read-outs to stabilise 
the detector (ranging from 10 to 34 frames depending on 
the expected source flux density) . The gain was 2 in most 
cases, but 1 in the case of sources that were expected to 
be relatively bright: HV12501 and 996, IRAS04496-6958, 
04545-7000, 05003-6712, 05112-6755, 05348-7024 and 
05568-6753, and WOH SG374. For most stars this re- 
sulted in a clear detection with S/N ratios of 10 to 100. In 
total 44 sources were observed with ISOCAM at 12 ^m. 

For sources that were expected to be stronger than 0.5 
Jy and which would therefore saturate the ISOCAM de- 
tectors with the LW10 filter, the 12 /zm photometry was 
obtained with the ISOPHOT instrument (Lemke et al. 
1996) using the 11.5 filter (~ IRAS 12 fim). These obser- 
vations were done using triangular chopping with a chop- 
per throw of 90". The aperture used for the observations 
was 52" in diameter. Integration times were 32 s on-source 
(and the same off-source), except for IRAS05294-7104 
that we integrated 64 s. A total of 13 sources were ob- 
served in this mode. For 53 sources we obtained PHOT-P 
photometry at 25 /im using the 25 filter (~ IRAS 25 /im) , 
triangular chopping with a chopper throw of 90", and an 
aperture of 52". Integration times ranged from 32 to 256 
s, depending on the expected flux density. In our Guar- 
anteed Time programme we finally observed 40 objects 
with ISOPHOT at 60 /im using the PHOT-C100 cam- 
era and filter 60 (~ IRAS 60 /im) and triangular chop- 
ping with a 150" chopper throw. Integration times ranged 
from 32 to 128 s, depending on the expected flux density. 
Unfortunately due to the reduced in-orbit sensitivity of 
the instrument and the problems with the calibration of 
the chopped measurements (especially for PHOT-C), we 
discovered after most observations had already been car- 
ried out that this was not the best observing strategy for 
the 60 /im photometry. Therefore, 7 objects were observed 
again in the Open Time using PHOT-C100 and filter 60 
in raster mapping mode, with 3x3 rasters and 45" raster 
steps in X and Y directions (spacecraft coordinates). The 
integration time per raster point was 128 seconds. 

In order to establish the carbon- or oxygen-rich na- 
ture of some of the programme stars we also obtained IR 
spectra for a number of them. In the Guaranteed Time 
15 objects were observed using PHOT-S in staring mode, 
with integration times of 256 or 512 s (1024 s for HV2379) 
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depending on the expected flux densities. The advantage 
of this instrument is that its spectral coverage is rather 
large (2 to 12 /im) at a reasonable resolution (~ 90). 
The sensitivity of the PHOT-S instrument, however, lim- 
its the detectability to sources with 12 /im flux densities 
above ~ 0.3 Jy. Furthermore, using staring observations 
the background cannot easily be determined. In this spec- 
tral region the diffuse emission is dominated by the zo- 
diacal emission, which, according to IRAS measurements, 
amounts to about <~ 0.1 Jy in the PHOT-S aperture at 10 

Considering this, we decided to obtain CAM-CVF 
spectra for 12 objects with a pixel field-of-view of 6" in 
beam LW-1. We did 25 exposures of 2.1 s each at gain 
2, after 50 read-outs to stabilise the detector. The un- 
precedented sensitivity of the ISOCAM instrument allows 
the observer to obtain spectra even for sources as faint 
as 100 mJy at 12 /im. Because of the long duration of a 
CVF observation, the spectral coverage chosen was only 
7 to 9.2 /im (with step 4) in LW-CVF1 and 9 to 14.1 
/zm (with step 2) in LW-CVF2, at a spectral resolution of 
~ 40. A big advantage of the CAM-CVF is that the spec- 
tra are obtained using an imaging technique. Therefore, a 
background spectrum was obtained simultaneously. These 
background spectra can be used to correct the PHOT-S 
spectra. We also obtained observations of 3 objects for 
which PHOT-S spectra had already been taken, in order 
to cross-check the results from the different instruments. 



4-1. Near-IR photometry 

Near-IR photometry was determined for each star at the 
time of the ISO observation, by interpolating ncar-IR 
lightcurves from our monitoring campaign at the South 
African Astronomical Observatory (SAAO) at Suther- 
land, South Africa (Whitelock et al., in preparation). 
Nearly always the lightcurve was sampled close in time 
to the ISO observation, but occasionally some extrapola- 
tion was necessary. The quoted uncertainties include an 
estimate, for each star, of the error introduced by the in- 
ter/extrapolation. For TRM45 and for the H-band magni- 
tude of IRAS05360— 6648 we have made use of the near-IR 
lightcurves and photometry presented by Wood (1998), 
after transformation to the SAAO system using Carter 
(1990). The near-IR photometry is listed in Table 4, along 
with the ISO photometry and the Julian Dates of the ISO 
spectroscopy. 

No near-IR counterparts could be identified with 
IRAS05568-6753 and 05289-6617. Two stars with near- 
IR colours much like those of unobscured M-type stars 
were monitored in the near-IR, but they show no variabil- 
ity. 




• = 12 fim (CAM) 
o = 12 /xm (PHOT) 
x = 25 iim (PHOT) 



■ 



2 

1.5 
1 

0.5 



0.2 0.4 0.6 O.B 1 



5 10 

F IS0 (Jy) 



Fig. 2. Comparison between ISO and our revised IRAS 
flux densities. 

5. ISO results and comparison with IRAS 
photometry 

The data were reduced using the PHOT and CAM Interac- 
tive Analysis software packages: PIA (Gabriel et al. 1997) 
version V7.1.2(c) and CIA (Ott et al. 1996) version V3.0, 
respectively. For a general description of the data and re- 
duction methods see the ISOPHOT Data Users Manual 
(Laureijs et al. 1998), and the ISOCAM Observer's Man- 
ual (1994) and ISOCAM Data Users Manual (Siebcnmor- 
gen et al. 1998). Details of the steps undertaken in re- 
ducing so-called Edited Raw Data (ERD) products to the 
finally derived flux densities and spectra can be found in 
Appendices A (photometry) and B (spectroscopy). The re- 
sulting ISO photometry is listed in Table 4, and the ISO 
spectra are presented in Figs. 4 & 5. 

The flux densities at 12 and 25 /im for the stars that 
were detected both by IRAS and ISO (Tables 3 & 4) are 
compared in Fig. 2. A bright regime where ISO and IRAS 
are consistent can be distinguished from a faint regime 
where ISO flux densities are systematically lower than 
IRAS flux densities. CAM is consistent with IRAS down 
to fainter levels (~ 0.2 Jy) than PHOT (~ 0.6 Jy). PHOT 
seems to undcr-estimate flux densities at levels between 
0.2 and 0.6 Jy by a factor -two. Below 0.2 Jy, both CAM 
and PHOT yield flux densities - 0.6 x IRAS. 

Flux density under-estimation may be caused by the 
difficulty of the detectors to respond to low signals. CAM 
12 /i flux densities below — 0.2 Jy may have been under- 
estimated as the K-a method does not adequately correct 
for non-stabilised signals if stabilisation is not reached well 
within less than half the integration time. The 12 and 25 
/im PHOT measurements, although employing different 
detector materials, show exactly the same trend. The 12 
(either CAM or PHOT) and 25 /im observations were per- 
formed in the same orbit. Default responsivities for the 
12 and 25 /im PHOT measurements were thought to be 
(much) lower in 1996 than in the current calibrations. 
Adopting those early values, the ISOPHOT photometry 
would be consistent with the IRAS photometry to a high 
degree. Subsequent revisions of the default responsivities 
have lead to higher values, approaching the values we ob- 
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tain from the (chopped) FCS measurements. The ratio of 
ISO and IRAS flux densities at 25 fiia is 0.60 ±0.06 for the 
24 sources with ISO measurements before orbit 190, and 
0.68 ± 0.09 for the 21 sources measured after orbit 190. 
These ratios are very similar, despite the large differences 
in median IRAS flux density between these two samples: 
0.41 and 0.14 Jy, respectively. 

However, the discrepancy between the ISO and IRAS 
data may not be as great as it appears if we take plausi- 
ble selection effects into account. The stars in our sample 
were largely selected on the basis of their IRAS flux den- 
sity, but many of them were only just detected by IRAS. 
It is therefore likely that they were near the maximum of 
their pulsation cycles at the time of the IRAS observa- 
tion. In contrast, they will have been at random phases 
when the ISO observations were made. This will lead to 
a systematic difference between the IRAS and ISO flux 
densities for faint sources. A similar effect may explain 
the discrepancy between the PHOT and CAM behaviour 
for sources with flux densities in the range 0.2 to 0.6 Jy, as 
the brighter sources were selected for measurement with 
PHOT and the fainter sources with CAM. Ground-based 
10 fim (N-band) magnitudes of a subset of our ISO targets 
were on average ~ 30% fainter than measured by IRAS at 
12 [im (Paper IV). Although we explained this in terms of 
differences between the N-band and IRAS 12 /urn filters, 
it may actually reflect the same discrepancy seen between 
the ISO and IRAS photometry. Variability cannot be the 
complete explanation, though: for instance, the sources 
IRAS04407-7000, 4516-6902 and 05003-6712 were all 
near the maxima in their K- and L-band lightcurves at 
the time of the ISO photometry, yet their PHOT 25 /im 
flux densities of 0.584, 0.380, and 0.210 Jy, respectively, 
are still fainter than the IRAS flux densities by factors of 
0.77, 0.69, and 0.64, respectively. Interestingly, Reid et al. 
(1990, their Figure 5) show that for Fi 2 , 2 5^0.3 Jy both the 
PSC and Schwering & Israel (1990) over-estimate flux den- 
sities for point sources in the LMC by typically 20 to 50%. 
They attribute this to source confusion resulting from the 
large beam width and crowdedness in typical LMC fields. 



Table 4. ISO 12, 25 and 60 /zm photometry (in Jy). The near-IR magnitudes for the ISO-epochs is deduced from light curves obtained at SAAO 
( JD — 2, 450, 000 = orbit + 38). Values in parentheses are l-cx errors. The last column indicates when a PHOT-S or CAM-CVF spectrum was taken. 



Star 




JD 


J[mag] 


H[mag] 


K[mag] 


L[mag] 


F 12 (CAM) 


Fi 2 (PHOT) 


F 25 (PHOT) 


F 6 o (chop) 


F( i0 (map) 


Spectrum 


GRV0519-6700 


318 


12.63(0.03) 


11.36(0.02) 


10.67(0.02) 




0.004(0.001) 












HV12070 




274 


10.03(0.02) 


8.98(0.02) 


8.68(0.02) 


8.24(0.03) 


0.043(0.004) 




-0.023(0.003) 












787 


10.30(0.10) 


9.20(0.10) 


8.80(0.10) 


8.30(0.10) 












CAM-CVF 


HV12501 




754 


8.80(0.05) 


8.00(0.05) 


7.75(0.05) 


7.30(0.05) 


0.185(0.005) 




0.067(0.018) 


0.260(0.108) 






HV2360 




754 


9.00(0.05) 


8.15(0.05) 


7.75(0.05) 


7.20(0.05) 


0.331(0.003) 




0.138(0.019) 


0.207(0.042) 






HV2379 




288 


12.68(0.03) 


12.05(0.04) 


11.40(0.04) 


9.70(0.10) 


0.038(0.002) 






-0.027(0.048) 




PHOT-S 






626 




13.20(0.40) 


11.90(0.20) 




0.035(0.001) 
















729 




13.20(0.40) 


11.90(0.20) 




0.032(0.001) 




0.014(0.005) 






CAM-CVF 


HV2446 




288 


10.25(0.02) 


9.23(0.02) 


8.80(0.02) 


8.35(0.04) 


0.066(0.001) 




0.043(0.005) 


0.109(0.056) 










622 


10.60(0.10) 


9.60(0.10) 


9.10(0.10) 


8.35(0.05) 










0.361(0.135) 








754 


10.10(0.02) 


9.04(0.02) 


8.67(0.02) 


8.06(0.02) 












CAM-CVF 


HV5870 




754 


9.34(0.02) 


8.40(0.02) 


8.10(0.02) 


7.60(0.04) 


0.269(0.002) 




0.092(0.008) 


0.090(0.355) 






HV888 




195 


8.10(0.01) 


7.20(0.01) 


6.89(0.01) 


6.46(0.02) 




0.713(0.037) 


0.201(0.012) 


-0.237(0.626) 




PHOT-S 


HV916 




209 


8.55(0.03) 


7.60(0.03) 


7.25(0.02) 


6.80(0.02) 


0.380(0.004) 




0.176(0.015) 


0.065(0.171) 






HV996 




163 


8.93(0.01) 


8.01(0.01) 


7.58(0.01) 


6.83(0.01) 


0.601(0.006) 




0.364(0.026) 


0.326(0.074) 




PHOT-S 


IRAS04286- 


-6937 


217 




13.00(0.05) 


11.25(0.05) 


9.10(0.05) 


0.136(0.002) 




0.059(0.013) 


-0.101(0.130) 






IRAS04374- 


-6831 


195 




13.60(0.05) 


11.40(0.04) 


8.80(0.03) 


0.185(0.003) 




0.086(0.027) 


0.228(0.136) 




PHOT-S 






701 




13.60(0.10) 


11.40(0.10) 


8.90(0.10) 










0.309(0.069) 




IRAS04407- 


-7000 


217 


10.18(0.03) 


8.92(0.03) 


8.18(0.03) 


7.30(0.05) 




0.962(0.024) 


0.584(0.031) 


0.125(0.088) 










605 


11.70(0.05) 


10.20(0.05) 


9.20(0.05) 


8.10(0.05) 










0.473(0.132) 




IRAS04496- 


-6958 


195 


13.00(0.05) 


10.90(0.05) 


9.50(0.04) 


7.70(0.04) 


0.269(0.002) 




0.126(0.010) 


0.252(0.154) 




PHOT-S 






605 


12.40(0.05) 


10.40(0.05) 


8.95(0.05) 


7.60(0.05) 












CAM-CVF 


IRAS04498- 


-6842 


217 


10.95(0.05) 


9.65(0.05) 


8.70(0.02) 


7.70(0.02) 




0.486(0.013) 


0.221(0.028) 


0.292(0.112) 






IRAS04509- 


-6922 


217 


11.50(0.05) 


10.00(0.04) 


9.15(0.04) 


8.10(0.04) 




0.303(0.011) 


0.202(0.028) 


-0.021(0.231) 






IRAS04516- 


-6902 


202 


10.50(0.03) 


8.96(0.03) 


8.24(0.03) 


7.25(0.03) 




0.854(0.026) 


0.380(0.051) 


0.349(0.135) 






IRAS04530- 


-6916 


202 


13.65(0.05) 


11.60(0.04) 


9.73(0.02) 


7.60(0.03) 




2.144(0.079) 


4.105(0.054) 


37.656(2.639) 






IRAS04539- 


-6821 


229 




14.30(0.05) 


11.80(0.04) 


8.80(0.10) 


0.244(0.002) 




0.089(0.013) 


0.106(0.053) 






IRAS04545- 


-7000 


195 




11.75(0.03) 


9.40(0.02) 


7.15(0.03) 


0.836(0.008) 




0.927(0.027) 


0.641(0.192) 




PHOT-S 






605 




12.80(0.10) 


10.30(0.10) 


8.10(0.10) 










0.286(0.206) 




IRAS04557- 


-6753 


229 




14.60(0.05) 


12.45(0.04) 


9.60(0.05) 


0.164(0.002) 




0.090(0.027) 


-0.072(0.137) 






IRAS05003- 


-6712 


229 


12.30(0.05) 


10.65(0.05) 


9.45(0.05) 


8.15(0.05) 


0.362(0.003) 




0.210(0.007) 


0.173(0.079) 




PHOT-S 


IRAS05009- 


-6616 


229 


14.70(0.10) 


12.65(0.05) 


10.70(0.05) 


8.45(0.05) 


0.284(0.003) 




0.082(0.017) 


-0.032(0.047) 






IRAS05112- 


-6755 


288 




14.70(0.10) 


12.00(0.05) 


8.80(0.05) 


0.414(0.004) 






0.078(0.133) 




PHOT-S 






629 




14.70(0.10) 


11.90(0.05) 


8.65(0.05) 


0.360(0.003) 




0.108(0.011) 






PHOT-S 


IRAS05113- 


-6739 


288 




14.35(0.05) 


12.12(0.03) 


9.06(0.02) 


0.316(0.003) 






0.398(0.115) 










629 




13.90(0.10) 


11.60(0.10) 


8.75(0.05) 


0.209(0.002) 
















729 




13.90(0.05) 


11.60(0.05) 


8.75(0.05) 


0.254(0.002) 




0.067(0.013) 








IRAS05128- 


-6455 


628 


13.60(0.10) 


12.10(0.10) 


10.55(0.05) 


8.55(0.05) 


0.226(0.002) 




0.061(0.008) 






PHOT-S 


IRAS05190- 


-6748 


288 






13.10(0.05) 


9.70(0.10) 


0.346(0.003) 




0.163(0.022) 


0.190(0.060) 




PHOT-S 


IRAS05289- 


-6617 


209 










0.157(0.002) 




0.202(0.009) 


0.382(0.049) 










668 


















0.281(0.069) 








729 




















CAM-CVF 



Table 4. (continued) The near-IR photometry of Wood (1998) is used for TRM45 and in part for IRAS05360— 6648 (H-band), after transformation to the 
SAAO system (Carter 1990). 



Star 




JD 


J[mag] 


H[mag] 


K[mag] 


L[mag] 


Fi 2 (CAM) 


IRAS05294- 


7104 


754 


11.80(0.05) 


10 


00 





05) 


8 


90 





05) 


7 


60 





05) 




IRAS05295- 


7121 


209 




13 


85 





05) 


11 


75 





05) 


9 


35 





05) 


0.143(0.001) 


IRAS05298- 


6957 


209 












11 


60 





20) 


8 


60 





20) 








729 


14.10(0.20) 


12 


50 





10) 


11 


10 





10) 


8 


50 





20) 




IRAS05300- 


6651 


209 




14 


70 





20) 


12 


20 





10) 


9 


10 





10) 


0.246(0.003) 


IRAS05329- 


6708 


163 


17.00(0.20) 


12 


70 





10) 


10 


40 





05) 


8 


25 





05) 




IRAS05348- 


7024 


202 




13 


60 





30) 


11 


60 





20) 


8 


70 





20) 


0.525(0.005) 






732 




13 


80 


'° 


30) 


12 


70 





20) 


9 


30 





20) 




IRAS05360- 


6648 


202 




15 


00 





40) 


12 


30 





05) 


9 


50 





10) 


0.171(0.002) 


IRAS05402- 


6956 


217 




13 


50 





10) 


10 


60 





05) 


8 


00 





04) 








662 


14.40(0.10) 


11 


45 





05) 


9 


40 





03) 


7 


15 





03) 








782 


14.50(0.10) 


11 


50 





05) 


9 


50 





05) 


7 


25 





05) 




IRAS05506- 


7053 


173 




16 


00 





30) 


13 


30 





10) 


10 


00 





10) 

J J 








228 




16 


40 





30) 


13 


60 





10) 


10 


30 





10) 




IRAS05558- 


7000 


209 


11.90(0.04) 


10 


10 





04) 


8 


90 





03) 


7 


60 





04) 








781 


13.50(0.05) 


11 


60 





03) 


10 


03 





01) 


8 


20 





05) 




IRAS05568- 


6753 


186 




























0.412(0.004) 


SHV0454030 


-675031 


318 


14.20(0.10) 


12 


67 





04) 


11 


85 





03) 










0.006(0.002) 


SHV0500193 


-681706 


229 


14.00(0.20) 


12 


00 





10) 


10 


70 





10) 


9 


50 





10) 


0.030(0.005) 






754 


13.10(0.10) 


11 


45 





05) 


10 


25 





05) 


9 


00 





05) 




SHV0500233 


-682914 


229 


13.05(0.05) 


11 


30 





05) 


10 


00 





05) 


8 


50 





05) 


0.076(0.001) 






729 


13.30(0.10) 


11 


40 





10) 


10 


20 





05) 


8 


60 





05) 




SHV0502469 


-692418 


304 


12.93(0.03) 


11 


61 





02) 


10 


81 





02) 


10 


30 





05) 


0.002(0.001) 


SHV0521050 


-690415 


788 


11.07(0.02) 


9 


77 





02) 


9 


23 





02) 


8 


40 





10) 


0.062(0.004) 


SHV0522023 


-701242 


788 


12.60(0.10) 


11 


60 


[0 


10) 


11 


40 





10) 










0.001(0.001) 


SHV0522118 


-702517 


217 


13.00(0.20) 


12 


00 





20) 


11 


00 





20) 


9 


70 





20) 


0.028(0.001) 


SHV0524565 


-694559 


217 


12.50(0.04) 


11 


40 





05) 


10 


75 





05) 










0.003(0.001) 


SHV0526001 


-701142 


217 


13.50(0.10) 


12 


00 





10) 


10 


45 





05) 


9 


15 





05) 


0.037(0.001) 


SHV0530323 


-702216 


788 


11.80(0.10) 


10 


80 





10) 


10 


35 





05) 










0.008(0.002) 


SHV0535442 


-702433 


782 


12.90(0.20) 


11 


50 





10) 


10 


65 


o 


05) 


10 


30 


o 


20) 


0.009(0.001) 


SP77 30-6 




195 


10.80(0.10) 


9 


60 





10) 


9 


20 





10) 


8 


50 





10) 


0.139(0.001) 






794 


10.30(0.10) 


9 


30 





10) 


8 


80 





10) 


8 


20 





10) 




TRM45 




788 


16.20(0.20) 


13 


60 





10) 


11 


55 





05) 


9 


50 





10) 


0.076(0.001) 


TRM72 




788 


15.22(0.20) 


12 


81 





03) 


10 


96 





02) 


8 


60 





05) 


0.161(0.001) 


TRM88 




788 


13.50(0.10) 


11 


78 





02) 


10 


24 





02) 


8 


55 





05) 


0.101(0.001) 


WBP14 




209 


13.35(0.05) 


11 


60 


;o 


05) 


10 


55 





03) 


9 


50 





10) 


0.022(0.001) 






616 


13.40(0.10) 


11 


60 


;o 


10) 


10 


55 





05) 


9 


50 





05) 




WOH G64 




229 


9.58(0.02) 


7 


96 





02) 


6 


98 





02) 


5 


32 





02) 




WOH SG374 


788 


9.90(0.02) 


9 


10 





02) 


8 


64 





02) 


7 


69 





05) 


0.490(0.005) 



Fi 2 (PHOT) F 25 (PHOT) F 60 (chop) F 60 (map) 



Spectrum 



0.680(0.045) 0.394(0.035) 
0.007(0.011) 
0.303(0.015) 0.359(0.014) 



0.442(0.018) 



0.455(0.020) 



-0.035(0.017) 
0.667(0.011) 



0.114(0.017) 
0.531(0.035) 
0.208(0.031) 

0.082(0.013) 
0.429(0.022) 



12.112(0.965) 



-0.011(0.014) 
0.003(0.021) 
0.007(0.002) 

-0.019(0.030) 
0.008(0.009) 
0.029(0.022) 
0.077(0.018) 

0.037(0.007) 
0.026(0.025) 
0.051(0.014) 
-0.005(0.011) 

13.786(0.354) 
0.188(0.020) 



0.079(0.097) 
0.233(0.074) 
0.414(0.279) 

0.068(0.049) 
0.479(0.082) 
0.170(0.108) 

0.174(0.067) 
0.944(0.: 



029(0 
517(0 


015) 
016) 


0.109(0.062) 
0.167(0.073) 
0.291(0.100) 


285(0 


008) 


0.383(0.073) 


030(0 


016) 


0.098(0.026) 


004(0 


033) 


-0.004(0.097) 



-0.245(0.527) 



0.170(0.125) 
0.051(0.063) 

0.225(0.073) 



PHOT-S 
CAM-CVF 

PHOT-S 

CAM-CVF 



CAM-CVf 

PHOT-S co 
CAM-CVP 

PHOT-S | 

a 
hi 

$ 

CAM-CVg 

CAM-CV^ 

O 
Cd 



CAM-CVF 



1.372(0.539) 

0.211(0.293) 
4.647(0.298) PHOT-S 
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Fig. 3. Histogram of the distribution of ISO 60 /im flux 
densities (chopped measurements). The dotted vertical 
line indicates 3-cr flux density derived from the distribu- 
tion of negative flux densities. All flux densities over 1 Jy 
are piled up in the last bin. 

A histogram of the distribution of ISO 60 /im flux den- 
sities (Fig. 3), leaving out the mapping observations, illus- 
trates the detection rate. Considering negative flux den- 
sities indicating non-detection, and assuming a Gaussian 
distribution around zero flux density for non-detections, 
we estimate a 1-cx detection to have 0.103 Jy. There are 
12 sources with flux densities exceeding 3-cr, i.e. proba- 
ble detections. This does not take into the account the 
large errors on some of the individual measurements, and 
a 3-cr detection may still turn out to be spurious (an ex- 
ample is WBP14). On the other hand, the distribution 
below 3-cr is certainly skewed towards positive flux densi- 
ties. Projecting the negative flux density distribution onto 
the positive domain, we estimate that there are probably 
17 more detections between and 3-cr, and a total of 14 
non-detections. 

In the IRAS 60 /im data we found 8 detections and 17 
tentative detections (Table 3). The 0.1 Jy assigned to the 
faintest IRAS 60 /im flux densities compares well with the 
ISO 1-cr detection threshold of the chopped measurements. 
Of the 8 IRAS 60 /im detections 6 have ISO chopped mea- 
surements, all of which yield higher flux densities than 
IRAS — by a factor 1.8 on average. This is in contrast 
to the 12 and 25 /im photometry, where ISO flux den- 
sities are generally lower than those measured by IRAS. 
The 14 ISO chopped measurements of IRAS tentative de- 
tections also yield higher flux densities than did IRAS - 
by a factor of 1.5 on average, although some individual 
ISO measurements are fainter than the IRAS ones. None 
of these ISO measurements is negative, indicating that 
many of the IRAS 60 tentative detections are indeed 
real. 

The 7 mapping observations all agree with the chopped 
measurements within 2-cr, although these errors can be 
large. There is no tendency for one of these two meth- 
ods to yield higher flux densities than the other. As we 
do not expect strong variability at 60 /im, which traces 
cool dust some distance from the stars, flux densities from 
mapping and chopped measurements are averaged. The 



error estimates of the mapping measurements are system- 
atically larger than those of the chopped measurements. 
This may be due to the fact that, for the mapping data, 
the flux density of the star was determined from the in- 
ner 3x3 pixels. The contribution of the background to 
these 9 pixels is considerable. Also, the reliability of the 
error estimate for the central pixel in the chopped data as 
produced by PIA is unknown. There is great difficulty in 
extracting reliable photometry and associated errors from 
either mapping or chopped measurements at 60 /im, for 
stellar sources in fields like the LMC. This is mainly due to 
the complex background and limited spatial resolution of 
PHOT at these wavelengths. IRAS05289-6617 has a very 
smooth background, being situated in the linc-of-sight to 
the supergiant shell LMC4 (Meaburn 1980). Indeed, ISO 
mapping and chopped measurements arc relatively pre- 
cise for this source, and agree nicely with the 60 (im flux 
density measured from the IRAS data. 

6. Discussion 

6.1. Chemical types from ISO spectroscopy 

The presence or absence in the ISO spectra (Figs. 4 & 5) of 
discrete dust emission and molecular absorption bands can 
be used to distinguish between carbon- and oxygen-rich 
circumstellar envelopes (e.g. Merrill & Stein 1976a, b,c). 
The results are summarised in Table 5. 

Amorphous oxygen-rich dust may give rise to strong 
and broad silicate emission between ~ 8 and 13 /im, peak- 
ing at ~ 9.7/im (the exact location may differ from this 
by a few tenths of /im). The late-M stars HV2446, 888, 
996, and SP77 30—6 have prominent silicate emission. In 
optically thick cases the silicate feature turns into absorp- 
tion. All spectra of OH maser sources show the silicate 
feature in self-absorption: IRAS04545-7000, 05298-6957, 
05329-6708, 05402-6956, and WOH G64. 

Oxygen-rich molecules do not provide clear diagnostics 
of the chemical type of CSEs at our signal-to-noise and 
spectral resolution. We already mentioned that shallow 
absorption around 3 /im in oxygen-rich sources is most 
likely due to an artifact in the responsivities, rather than 
H 2 ice. 

Crystalline carbon-rich dust sometimes gives rise to a 
SiC (graphite) emission feature peaking at ~ 11.3/zm, and 
narrower than the silicate feature. The CVF spectrum of 
IRAS05289-6617 (Fig. 4) shows the best example of this. 

Carbon-rich molecules have several strong absorption 
bands in our spectral region, all from HCN and C2H2. The 
strongest is at 3.1 /im, but the problem with the respon- 
sivities limits the number of unambiguous detections to 
one (IRAS04496— 6958). Related, but weaker, absorption 
is visible at 3.8 /im. More absorption bands are located 
around 5, 8 and 14 /im. Unfortunately, the 5 /im band 
falls entirely in the blind spectral region of PHOT-S. The 
8 and 14 /im bands are at the edges of the CVF spectra 
and hence difficult to identify. 
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Fig. 5. The PHOT-S spectra of obscured AGB stars (and RSGs) in the LMC. Emission and/or absorption centred at 
~ 9.7 yum suggests oxygen-rich dust (e.g. HV888 and IRAS04545— 7000). Absorption at 3 is seen in carbon star 
photospheres (e.g. IRAS04496— 6958), but artifacts in the PHOT-S rcsponsivities also mimic weak depression at 3 /j,m 
in the spectra of unambiguous oxygen-rich stars (e.g. IRAS04545— 7000). 



Some other spectra show merely a featureless dust con- filters are unknown, we adopt here the IRAS zero-points, 
tinuum around 10 /im. Best examples are the CVF spec- This results in the following definitions for the (not colour- 
trum of SHV0500193-681706 and the PHOT-S spectrum corrected) mid-IR magnitudes 
of IRAS05568— 6753. These spectra suggest pure amor- 
phous carbon dust emission. I 12 ] = -2.5 log(Fi 2 /28.3) (1) 



6.2. IR colour-colour diagrams 

The ISO 12, 25 and 60 /xm filters are similar but not iden- 
tical to the IRAS filters. As the zero-points of these ISO 



[25] = -2.51og(F 25 /6.73) 
[60] = -2.51og(F 60 /1.19) 



(2) 
(3) 
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Fig. 5. (continued) 



6.2.1. Diagram of (K - [12]) versus (H - K) 

The (K — [12]) versus (H — K) colour-colour diagram sep- 
arates carbon- from oxygen-rich stars in samples of ob- 
scured stars in the MCs (Papers II & IV). Indeed, the 
distributions of carbon- and oxygen-rich stars using ISO 
and SAAO photometry define clear sequences in this di- 
agram (Fig. 6). The sequences are fit by eye, with the 
carbon sequence the same as in Paper IV: 



(H — K) = 0.36 x(K- [12]) 



(4) 



but the oxygen sequence a simple, yet somewhat steeper 
function than in Paper IV: 



(H — K) = 0.3 + 0.0003 x (K - [12]) E 



(5) 



Although the stars with spectral type M follow the oxy- 
gen sequence very well, the carbon stars show a large scat- 
ter around the carbon sequence with several carbon stars 
on or beyond the region populated by M stars, at small 
(H — K) but large (K — [12]) magnitudes. This scatter 
contrasts with the tight carbon sequence that is observed 
in the Milky Way (Fig. 3 in Paper IV). We suspect that 
this is in part caused by the severe crowding in some LMC 
fields, affecting the near-IR aperture photometry. Differ- 
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Table 5. Chemical types. Optical spectra (Opt Sp) in- 
clude objective prism and CCD spectroscopy up to ~ 
1/xm. ISO spectroscopy (ISO Sp) comprises PHOT-S and 
CAM-CVF observations. IR colour-colour diagrams (IR 
col) can in some cases be reasonably conclusive too: we 
here use (K - [12]) and ([12] - [25]) versus (K - L) dia- 
grams. At radio wavelengths, OH, SiO and/or H2O maser 
emission is detected from some oxygen-rich sources. 



0.06 
3 0.04 
6^0.02 

0.1 

0.05 
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GRV0519-6700 C 
HV12070 MS3/9 oxygen? 

HV12501 
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Fig. 4. The CAM-CVF spectra of obscured AGB stars 
in the LMC. Open symbols represent spectro-photometric 
points that are prone to have flux densities that are over- 
estimated due to stabilisation difficulties. The spectral 
shape is best represented by the solid symbols (squares 
for the short-, disks for the long- wavelength region). Emis- 
sion and/or absorption centred at ~ 9.7 fxm is indica- 
tive of oxygen-rich dust (e.g. IRAS05402-6956 and SP77 
30—6), whereas carbon-rich dust may show emission at 
- 11.3 [im (e.g. IRAS05289-6617). A featureless contin- 
uum around 10 /xm also strongly suggests carbon-rich dust 
(e.g. SHV0500193-681706). 

Table 5. (continued). 
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WOH SG374 
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ences in the strength of absorption in the H-band by car- 
bonaceous molecules may cause additional scatter among 
carbon stars (Bessell & Wood 1983; Catchpole & White- 
lock 1985). 
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Fig. 6. (K — [12]) versus (H — K) diagram. Stars are dis- 
tinguished by their chemical types inferred from spectro- 
scopic methods: carbon stars (solid disks), M stars (open 
disks), MS stars (open squares), S stars (open triangles), 
and stars of which the chemical type is a priori unknown 
(crosses) . Carbon stars and oxygen stars define sequences 
in this diagram, indicated by a dotted and solid curve, 
respectively. 

6.2.2. Diagram of (K - [12]) versus (K - L) 

The (K — [12]) versus (K — L) colour-colour diagram shows 
much less scatter around well-defined carbon and oxygen 
sequences (Fig. 7). This makes it a much more powerful 
diagnostic diagram than the (K — [12]) versus (H — K) 
diagram in typifying the chemical composition of the cir- 
cumstellar dust. Noguchi et al. (1991a) introduced a very 
similar diagnostic using (L — [12]) and (K — L) colours. 
We note, however, that some of the peculiar stars in our 
(K— [12]) versus (H—K) diagram were too blue and hence 
too faint to be detected in the L-band. Still, the tight se- 
quences prove that both the SAAO and ISO photometry 
are of good quality when comparing individual stars. We 
fit (by eye) a linear carbon sequence: 



(K-L) = ^x(K-[12])-± 



(6) 



and a superposition of even polynomials for the oxygen 
sequence: 

(K—L) = 0.35+0.007x(if-[12]) 2 +0.0014x(if-[12]) 4 (7) 

6.2.3. Diagram of ([12] - [25]) versus (K - L) 

Another colour-colour diagram that separates carbon- 
from oxygen-rich sources is the ([12] — [25]) versus (K — L) 
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Fig. 6. 



(K—L) 

(K — [12]) versus (K — L) diagram. Symbols as in 



diagram (Fig. 8). The confirmed oxygen- rich sources show 
a linear relationship between the ([12] — [25]) and (K — L) 
colours, possibly flattening out at (K — L) > 1.5 mag. The 
LMC stars generally follow the separation determined for 
galactic stars (dotted line in Fig. 8, taken from Epchtein 
et al. 1987). 

6.2.4. Diagram of ([25] - [60]) versus ([12] - [25]) 

The aim of obtaining 60 /im flux densities for stars in the 
LMC is mainly to probe the coolest circumstellar dust. 
The 60 /zm flux density is expected to increase as pro- 
longed mass loss first extends the CSE and again as re- 
duced mass loss results in a detached shell. This evolution 
might be seen in ([25] — [60]) versus ([12] — [25]) diagrams 
(Fig. 9, see also van der Veen & Habing 1988). Unfor- 
tunately, the accuracy of the ISO photometry at 60 /im 
is not very high for most of these LMC sources, and the 
diagram contains a lot of scatter. 

Perhaps the most obvious thing to learn from this di- 
agram is that carbon stars tend to be relatively bright 
at 60 (jm, yielding ([25] - [60]) - 1.5 to 3 mag. Al- 
though oxygen-rich sources can have similar colours, there 
are many oxygen-rich sources with ([25] — [60]) < 2 and 
([12] — [25]) > 0.6 mag, colours not seen for any carbon 
star in our sample. This is similar to the findings of van der 
Veen & Habing (1988), but our LMC sources have bluer 
([12] - [25]) and redder ([25] - [60]) colours than do their 
Milky Way sources. However, the LMC ([12] - [25]) colours 
do not differ much from those discussed by Le Bertre et 
al. (1994). 
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Fig. 8. ([12] — [25]) versus (K — L) diagram. Symbols are 
as in Fig. 6. Oxygen-rich sources and carbon stars occupy 
distinct areas in this diagram. The dividing line (dotted) 
between stars with carbon- and oxygen-rich dust is taken 
from Epchtcin ct al. 1987. 

6.3. Comments on particular objects 

6.3.1. GRV0519-6700 

The referee Dr. Peter Wood conveys that an optical spec- 
trum of GRV0519— 6700 shows it to be a carbon star, in 
good agreement with its IR colours of (H — K) = 0.7 and 
(K - [12]) = 1.0 mag. 

6.3.2. HV12070 

The CVF spectrum of HV12070 shows only a hint of the 
silicate feature, whilst the IR colours cannot distinguish 
between oxygen- and carbon-rich dust of the optically thin 
CSE of this MS-type star. 

6.3.3. HV2379 

The PHOT-S spectrum of HV2379 suggests SiC emission, 
but its CVF spectrum does not. This may be a result of 
changes in the properties of the CSE or the dust. Its IR 
colours leave no doubt about the carbon- rich nature of the 
dust. 

6.3.4. HV2446, 888, 996, and SP77 30-6 

These late-M stars all have prominent silicate emission 
and IR colours that unambiguously indicate oxygen-rich 
dust. 
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Fig. 9. ([25] - [60]) versus ([12] - [25]) diagram. Symbols 
are as in Fig. 6. Carbon stars are not well separated from 
oxygen-rich sources, although carbon stars seem to be rel- 
atively bright at 60 /um. 



6.3.5. IRAS04286-6937, 04539-6821, 04557-6753, 
05009-6616, 05113-6739, 05295-7121, 
05300-6651, 05360-6648, and TRM45 and 72 

The position of these objects in the (K — [12]) versus 
(H — K) or (K — L) colour-colour diagrams does not clarify 
the chemical composition of their CSEs. The ([12] — [25]) 
versus (K — L) diagram, however, unambiguously indi- 
cates that the dust around these stars is carbon rich. The 
IR colours of IRAS05113-6739 at the three ISO epochs for 
this star all lie along the carbon sequences in the (K — [12]) 
versus (H — K) and (K — L) diagrams. Ground-based L- 
band spectra of IRAS05009-6616 and 05300-6651 show 
the 3.1 ^im absorption feature due to HCN and C2H2 
molecules, indicating carbon-rich photospheres (van Loon 
et al. 1999). 

6.3.6. IRAS04374-6831 

The position of IRAS04374-6831 in the ([12] -[25]) versus 
(K — L) diagram indicates carbon-rich dust. Its PHOT-S 
spectrum, which does not clearly reveal the chemical com- 
position of the dust by itself, is then marginally consistent 
with SiC emission. 

6.3.7. IRAS04496-6958 

IRAS04496— 6958 shows strong absorption by carbona- 
ceous molecules at 3.1 fim, already known from ground- 
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based L-band spectroscopy (van Loon et al. 1999). Re- 
lated, but weaker, absorption is visible at 3.8 /<m, and 
possibly around 8 /zm. Surprisingly, this carbon star has 
silicate emission too, indicating the presence of oxygen- 
rich dust (see Trams et al. 1999). Its IR colours indicate 
carbon- rich dust, hence the oxygen-rich dust is only a mi- 
nor component. 

6.3.8. IRAS04530-6916 

With (K-L) = 2.13, {K - [12]) = 6.9 and ([12] - [25]) = 
2.3 mag, the IR colours of IRAS04530-6916 imply that 
the dust around this very luminous and red object must 
be oxygen rich. 

6.3.9. IRAS04545-7000, 05298-6957, 05329-6708, 
05402-6956, and WOH G64 

These OH maser sources all show the silicate feature in 
self-absorption, and also their IR colours clearly indicate 
oxygen-rich dust. 

6.3.10. IRAS05003-6712 

The IR colours of IRAS05003— 6712 unambiguously clas- 
sify the dust as oxygen rich. The PHOT-S spectrum shows 
a hint of the silicate feature. A ground-based L-band spec- 
trum of this star shows a featureless continuum around 3.1 
fim, indicating an oxygen-rich photosphere (van Loon et 
al. 1999). 

6.3.11. IRAS05112-6755 

The dust around IRAS05112— 6755 is classified as carbon 
rich on the basis of the position in the ([12] — [25]) versus 
(K — L) diagram. There is a hint of 8 /im absorption in the 
PHOT-S spectrum of IRAS05112-6755. A ground-based 
L-band spectrum of this object shows the strong absorp- 
tion at 3.1 fim found in carbon-rich stellar photospheres 
(van Loon et al. 1999). 

6.3.12. IRAS05128-6455 and 05190-6748 

The absence of clear indications for the presence of the sil- 
icate feature in the PHOT-S spectra of these stars suggest 
that their dust may be carbon rich, which is also indicated 
by their ([12] - [25]) and (K - L) colours. 

6.3.13. IRAS05289-6617 

The CVF spectrum of IRAS05289-6617 shows prominent 
SiC emission. Hence it is probably a mass-losing carbon- 
rich AGB star in the LMC rather than a foreground ob- 
ject. We have not yet identified its near-IR counterpart. 



6.3.14. IRAS05348-7024 

The CVF spectrum of IRAS05348-7024 shows weak SiC 
emission. The carbon-rich nature of the dust around this 
object is also indicated by its position in the ([12] — [25]) 
versus (K — L) diagram. 

6.3.15. IRAS05506-7053 

IRAS05506— 7053 is the only star in our sample that could 
not be detected at 12 ^im. Assuming a 12 /im flux density 
< 0.03 Jy, the (K — [12]) colour would be < 6.2 mag and 
probably ([12] - [25]) > 1.5 mag. At (K - L) = 3.3 mag, 
this suggests an oxygen-rich CSE. 

6.3.16. IRAS05558-7000 

The CVF spectrum of IRAS05558-7000 is similar to the 
CVF spectra of IRAS05298-6957 and 05402-6956, show- 
ing silicate emission that is becoming optically thick at 10 
^m. The IR colours of IRAS05558-7000 unambiguously 
imply that the dust is oxygen-rich. 

6.3.17. IRAS05568-6753 

The PHOT-S spectrum of IRAS05568-6753 shows a fea- 
tureless dust continuum around 10 /xm, suggesting pure 
amorphous carbon dust emission. The near-IR counter- 
part of this object has yet to be identified. 

6.3.18. SHV0500193-681706 

The CVF spectrum of SHV0500193-681706 shows a fea- 
tureless dust continuum around 10 /im, suggesting pure 
amorphous carbon dust emission. The carbon-rich nature 
of the dust is confirmed by the position in the (K — [12]) 
versus (K — L) diagram. Inaccuracy of its 25 ^.m flux den- 
sity causes the rather odd position among the oxygen-rich 
stars in the ([12] — [25]) versus (K — L) diagram. 

6.3.19. SHV0500233-682914 

The CVF spectrum of SHV0500233-682914 shows a hint 
of SiC emission, and also its IR colours clearly indicate 
that the dust around this star is carbon rich. 

6.3.20. SHV0502469-692418, 0522023-701242 and 
0524565-694559 

The carbon star SHV0502469-692418, the M-type 
star SHV0522023-701242 and the MS-type star 
SHV0524565— 694559 are surrounded by an optically 
thin CSE and hence it is difficult to classify the chemical 
type of their dust from IR colour-colour diagrams. 
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6.3.21. SHV0522118-702517 

SHV0522118-702517 was tentatively classified an S-type 
star by Hughes & Wood (1990). Its IR colours are clearly 
similar to those of carbon stars. This suggests that carbon- 
rich dust dominates the absorption and emission charac- 
teristics of the CSE despite the under-abundance of car- 
bon atoms in its photosphere. Noguchi et al. (1991b) show 
that the IR colours of the CSE indicate oxygen-rich dust 
in case of an MS-type star. Also, CS stars show 3 fim 
absorption from HCN and C2H2 molecules, whereas SC 
stars do not (Catchpole & Whitelock 1985; Noguchi & 
Akiba 1986). This suggests that carbon chemistry is dom- 
inant in CS stars, but not in SC stars. Thus, we identify 
SHV0522118-702517 with a CS star. Dust-enshrouded S 
stars — including MS and CS stars — that have {K — L) > 
1 mag are very rare in the Milky Way, and none are known 
with {K — L) > 2 mag (Noguchi et al. 1991b). Hence, with 
{K — L) = 1.3 mag, SHV0522118-702517 is among the 
most obscured S stars known. 

6.3.22. SHV0530323-702216 

The late-M type star SHV0535442- 702433 has ([12] - 
[25]) = 1.56 mag. Its near-IR colours are rather blue and 
{K — L) is not expected to be larger than unity. Hence the 
position of this object in the ([12] — [25]) versus {K — L) 
diagram suggests that the dust is oxygen rich. 

6.3.23. SHV0535442-702433 

The carbon star SHV0535442-702433 is surrounded by an 
optically thin CSE, and hence the IR colours are difficult 
to use for classifying the chemical type of the dust. The 
location among oxygen-rich stars in the ([12] — [25]) versus 
{K — L) diagram is caused entirely by the inaccuracy of its 
25 ^m flux density yielding a spuriously red ([12] — [25]) ~ 
3 mag. 

7. Conclusions 

ISO spectroscopy is used to determine the chemical type of 
the dust around obscured cool evolved stars in the LMC. 
ISO photometry at 12, 25 and 60 /im is presented, to- 
gether with quasi-simultaneous near-IR photometry from 
the ground (SAAO). The accuracy and sensitivity of the 
ISOPHOT photometry is not much better than can be 
achieved from properly treated IRAS data. The ISOCAM 
photometry is much more reliable because it is based 
on imaging, and an order of magnitude more sensitive 
than was IRAS. Colour-colour diagrams prove that rel- 
ative photometry is reliable. A combination of {K — [12]) 
and ([12] — [25]) versus {K — L) diagrams provide a reli- 
able way of distinguishing between carbon- and oxygen- 
rich dust, provided the CSE has sufficient optical depth. 
The combination of ISO spectra and photometry enabled 
us to securely classify the chemical type of the dust around 



nearly all stars in our sample. This was previously known 
for only a minority of the stars. Surprisingly, the {K — [12]) 
versus {H—K) diagnostic diagram contains a lot of scatter 
especially among carbon stars. 

Many of the obscured AGB stars in our sample arc 
carbon stars: 46% amongst the LMC stars that were de- 
tected by IRAS (Table 1). M stars were always found to 
be surrounded by oxygen-rich dust. In particular, all de- 
tected OH maser sources show self-absorbed silicate emis- 
sion. As in the Milky Way, the fact that no M star with 
carbon-rich dust has ever been found suggests that HBB 
cannot efficiently turn carbon stars back into oxygen- 
rich stars. The dust around the dust-enshrouded S star 
SHV05221 18-702517 has the characteristics of carbon- 
rich material, suggesting it is actually a CS star. Surpris- 
ingly, the dust around the carbon star IRAS04496— 6958 
has a (minor) oxygen- rich component (Trams et al. 1999). 
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Appendix A: ISO photometry 

A.l. CAM 12 fim imaging-photometry 

The ISOCAM images were corrected for the dark image valid 
for the corresponding revolution and orbital position of the 
spacecraft, and corrected for glitches using the multiresolution 
median transform method. Although the current dark subtrac- 
tion algorithm produces satisfactory results, we nonetheless 
also applied our own IDL routine that optimises the dark sub- 
traction by appropriately scaling the dark and flatfield images 
in a rectangular annulus leaving a 12 x 12 pixel area centred 
at the stellar position on the array. The dark subtraction was 
improved slightly in some cases, whereas for the majority of 
the observations there was no significant difference. For each 
pixel, we applied a n-a rejection criterion to the evolution of 
the signal in time: values that deviate more than 2 a from the 
time-average were rejected. The images were divided by the de- 
fault calibration flatfield image to correct for the array pattern 
in the pixel responsivities, and time-averaged. 

Point source photometry was performed on the final image, 
using our own procedures written within the Munich Interac- 
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tive Data Analysis Software (MIDAS). We measured the flux 
density within a circular area (software aperture) around the 
stellar position, and subtracted the background flux density 
level determined from a concentric circular annulus between 
12 and 14 pixels radius from the stellar position — exclud- 
ing the pixels near the edge of the 32 x 32 pixels array where 
vignetting is evident. We repeated this for different software 
aperture radii, to create a magnitude profile (MP) of the stel- 
lar flux density versus radius. We used the 23 brightest stars 
— with flux densities ranging from 0.06 to 0.8 Jy — to create 
a template MP. Comparison of the individual MPs of the stars 
with the template yielded differential MPs. Where the differ- 
ential MP is constant with radius, reliable differential magni- 
tudes can be determined. We adopt the standard deviation in 
these points as 1-a errors on the differential photometry. The 
strength of the MP-method is the estimation of photometric 
errors, as well as the selection of the part of the MP that best 
resembles the PSF. Calibration of the template MP was done 
using the calibration conversion 4.13 ADU s^ 1 mJy- 1 (OLP V7 
calibration) . The template MP was consistent with a synthetic 
MP created from the known PSF, and indicated an uncertainty 
of ~ 1 mjy. 

All target stars were detected unambiguously, except for 
SHV0522023- 70124 which was marginally detected. Most were 
well centred — within about 10" — but IRAS05113-6739 
(last epoch only), IRAS05295-7121 and TRM45 were off 
by 20 to 30". Seven stars were accompanied by generally 
weaker additional point sources within the image: HV12501, 
SHV0500193-681706, 0500233-682914, 0502469-692418, and 
0535442-702433, and TRM72 and 88. All of these stars are 
either HV, SHV, or TRM sources, and the nature of the addi- 
tional point sources remains a mystery. 

A. 2. PHOT-P 12 and 25 [im chopped measurements 

The measurements each consist of 4 Destructive Readouts 
(DRs), with 127 Non-Destructive Readouts (NDRs) per ramp. 
All DRs and the first 8 NDRs per ramp were discarded, because 
these read-outs are considered unreliable due to the strong im- 
pact of the detector reset. The signal corresponding to the 
source intensity is derived from the slope of the integration 
ramp. The ramps were linearised to correct for non-linearities 
in the CRE (Cold-Read-Out Electronics) output voltages and 
debiasing effects, which mainly affect the long wavelength de- 
tectors. We attempted to remove glitches resulting from the im- 
pact of high energy particles by applying the two-threshold me- 
dian filtering technique, using distributions of 32 data points, 
three iterations, and thresholds of three standard deviations 
for flagging and re-accepting. The ramps were subdivided into 
four sections, so as to be more selective in discarding erroneous 
read-outs. The first 50% of the signals per chopper plateau 
were discarded to enable the detector response to stabilise at 
the level of the source flux density. The chopper plateaus were 
then treated separately in deglitching the signals. This was 
done by running a bin with a width of 8 signals over the chop- 
per plateau, taking individual steps, and iterating twice: after 
maximum/minimum clipping, signals were discarded if they 
were off by more than three standard deviations. A correc- 
tion was applied to account for the signal dependence on the 
reset interval which was used for the read-out sampling. The 



expected dark current for the orbital position of the space- 
craft was subtracted from the data. The differences between 
on-source and interpolated background signals were corrected 
for signal losses due to the rapid chopping. 

The main source of concern is that the detector response 
is not constant over time, and depends on the history of the 
detector illumination and on the levels of source and back- 
ground. Therefore, careful modeling of the detector behaviour 
is needed to determine the true source-background signal. At 
present the quality of our data and the degree of understanding 
of the detector characteristics is insufficient to allow for such 
an advanced data reduction. We opted instead for a simpler ap- 
proach, that also conforms better with the way the standard 
star observations and flux-density-calibrations were done. 

The in-band power was calibrated using the responsivity 
as derived from an FCS1 internal calibrator measurement re- 
duced in the same way as the scientific measurement, except 
that background subtraction and chopper frequency correction 
are not applicable. The ratio of the responsivities determined 
from FCS1 and default varies between 0.8 and 2.8, averaging 
1.4 around orbit 200 and 2.2 around orbit 700. We therefore 
feel that the FCS1 values reveal the complex detector- history 
dependence of the responsivity, and hence should be applied to 
calibrate the in-band power, instead of using the recommended 
default responsivity for chopped observations. 

Finally the median flux density was extracted, correcting 
for the point source flux density outside of the aperture. The 1- 
a errors are estimated by quadratic summation of the flux den- 
sity uncertainty and the error in the responsivity value. The 25 
^m flux densities of IRAS05113-6739, IRAS05295-7121 and 
TRM45 have almost certainly been under-estimated due to off- 
centering of the source with respect to the ISOPHOT aperture 
(see the previous subsection about the CAM imaging). 

A. 3. PHOT-C 60 [im chopped measurements 

Observations at a wavelength of 60 /im are severely hampered 
by the complex and bright sky background compared to the 
flux density of a typical point source in our sample. The sky 
in the direction of the LMC varies on a spatial scale of a few 
arcminutes, due to the presence of molecular cloud complexes 
within the LMC which are much brighter than circumstellar 
envelopes at 60 /im. This scale is not much larger than the 
PSF for a telescope under 1 m diameter. 

The measurements each consist of 4 DRs (8 for WOH G64) , 
with the number of NDRs per ramp depending on the expected 
flux density. Usually we chose 15, 31 or 63 NDRs per ramp, but 
for IRAS04530-6916 and 05329-6708 and for WOH G64 only 
7 and 3, respectively, were chosen. All DRs were discarded, as 
well as the first 3 or 6 NDRs per ramp in case of 31 or 63 NDRs 
per ramp, respectively (otherwise only one NDR per ramp was 
discarded). The ramps were linearised. If a sufficient number 
of data points were available, glitches were removed by apply- 
ing the two-threshold median filtering technique, using distri- 
butions of 16 data points, three iterations, and thresholds of 
three standard deviations for flagging and re-accepting. The 
ramps were subdivided into four sections, provided there were 
at least four data points per subramp. The first 50% of the 
signals per chopper plateau were discarded. The signals were 
deglitched by running a bin with a width of 16 signals over the 
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chopper plateau, taking individual steps, and iterating twice: 
after maximum/minimum clipping, signals were discarded if 
they were off by more than three standard deviations. A reset- 
interval correction was applied. The expected dark current for 
the orbital position of the spacecraft was subtracted from the 
data. A correction for vignetting of the 3x3 pixels C100 array 
was applied. The differences between on-source and interpo- 
lated background signals were corrected for signal losses due 
to the rapid chopping. This correction is very large — a factor 
of 2.75 — and rather uncertain (35%). 

The in-band power was calibrated using the responsivity 
as derived from an FCS1 internal calibrator measurement done 
immediately after the scientific measurement. The FCS respon- 
sivities were found to resemble the pixel response pattern in 
the scientific measurement better than the default responsivi- 
ties did. The FCS measurement was reduced in the same way as 
the scientific measurement, except that vignetting correction, 
background subtraction and hence a chopper frequency cor- 
rection are not applicable. The ratio of the FCS1 and default 
responsivities varied between 0.8 and 2.0, around a median of 
1.3 and with a standard deviation of 0.3. It varied from orbit 
to orbit and within an orbit, but with no clear time depen- 
dence — except that it was larger after orbit 190 than before. 
This behaviour closely resembles that seen for the PHOT-P 
detectors. 

Finally the median flux densities of the pixels were ex- 
tracted. The median value of the eight pixels surrounding the 
central one was subtracted from the central pixel value. Know- 
ing that the central pixel contains 66.35% of the total flux 
density of a well-centred point source, and the entire C100 ar- 
ray contains 91.75%, a correction factor of 1.5835 was applied. 
The (internal) 1-a errors are estimated by quadratic summa- 
tion of the error in the value for the central pixel and 1/V8 
times the median of the differences between background pixels 
and their median. This does not take into account the uncer- 
tainty in the absolute flux-density calibration. The flux density 
of IRAS05295— 7121 has almost certainly been under-estimated 
due to off-centering of the source with respect to the ISOPHOT 
aperture (see the earlier subsection about the CAM imaging) . 

A. 4- PHOT-C 60 \im mapping observations 

The measurements each consist of 64 DRs, with 63 NDRs per 
ramp (128 DRs and 31 NDRs per ramp for IRAS05402-6956 
and WBP14). All DRs were discarded, as well as the first 3 
or 6 NDRs per ramp in case of 31 or 63 NDRs per ramp, re- 
spectively. The ramps were linearised. Glitches were removed 
by applying the two-threshold median filtering technique, using 
distributions of 32 data points, three iterations, and thresholds 
of three standard deviations for flagging and re-accepting. The 
ramps were subdivided into four sections. The first 4 read-outs 
per raster point were discarded. The signals were deglitched 
by running a bin with a width of 64 signals through the raster 
point, taking individual steps, and iterating twice: after max- 
imum/minimum clipping, signals were discarded when they 
were off by more than three standard deviations. We applied 
the stability recognition method to derive the signal most likely 
to be near the true, stabilised one: a bin with a width of 64 
signals was run over the measurement, with intervals of 32 sig- 
nals, and only data within a confidence level of 0.95 were kept. 



A reset-interval correction was applied. The expected dark cur- 
rent for the orbital position of the spacecraft was subtracted 
from the data. The 3x3 pixels array was corrected for vi- 
gnetting. 

The in-band power was calibrated by interpolating the re- 
sponsivities derived from FCS1 internal calibrator measure- 
ments made before and after the scientific measurement. Use 
of the two FCS measurements was found to correct (in part) for 
drift behaviour of the signal, which was not the case if default 
responses were used. The FCS measurements were reduced in 
the same way as the scientific measurement, except that a vi- 
gnetting correction is not applicable and 16 data points were 
used for signal deglitching and stability recognition. Finally 
the median flux densities of the pixels were extracted. 

The values of the eight pixels surrounding the central one 
were replaced by their median value and added to the central 
pixel value. The values of the sixteen pixels along the rim of 
the map were also replaced by their median value, and this was 
taken to be the sky background to be used for correcting the 
total flux density in the inner 3x3 part of the map. Knowing 
that the central pixel contains 66.35% of the total flux density 
of a well-centred point source, and the entire C100 array con- 
tains 91.75%, the outer sixteen pixels contain a fraction of the 
flux density that is somewhere between and 8.25%. Assum- 
ing this fraction to be 4.125%, a factor of 1.1182 is derived to 
correct the background subtracted flux density for the stellar 
flux density contained in the outer sixteen pixels. The (inter- 
nal) 1-a errors are estimated by quadratic summation of the 
error in the value for the central pixel, 8/v8 times the me- 
dian of the differences between the eight surrounding pixels 
and their median, and 9/\/l6 times the median of the differ- 
ences between the outer sixteen pixels and their median. The 
error in the value for the central pixel was determined by the 
median of the differences between the values for the pixels that 
were centred on the star over the course of the mapping, and 
their median. The error estimate does not take into account 
the uncertainty in the absolute flux-density calibration. 

Appendix B: ISO spectroscopy 

B.l. CAM-CVF spectro-photometry 

The CAM-CVF cube (X, Y, A) was corrected for the dark im- 
age using a model for the dependence on the revolution and 
orbital position of the spacecraft, and corrected for glitches by 
applying a multiresolution median transform. For each pixel, 
we applied a K-a rejection criterion to the evolution of the sig- 
nal in time: values that deviate more than 2 a from the time- 
average were rejected. Attempts to correct for the transient 
behaviour of the signal by applying various models yielded re- 
sults no better than this method. The images were divided by 
the default calibration Hatfield image to correct for the array 
pattern in the pixel responsivities, and time-averaged. 

The spectra were constructed by obtaining photometry 
from the images that each correspond to a different position 
of the CVF. For each target we integrated the signal over the 
pixels that were significantly above the background level. The 
number of pixels was limited to avoid sampling excessive back- 
ground, and ranged from 1, for the faintest, to 9, for the bright- 
est sources. The short- and long- wavelength parts of the CVF 
were treated separately because the star was often centred at 



20 



Trams et al.: ISO observations of AGB stars in LMC 



a slightly different position on the array. The background level 
was determined by taking the median value within a three- 
pixel wide circular annulus around the star, with an inner ra- 
dius of three pixels. The stellar flux density was corrected for 
the wavelength dependence of the PSF. This correction was 
derived by doing photometry on the PSF with equivalent in- 
tegration areas but inversely proportional to the wavelength. 
The first few (~ 4) spectro-photometric points are prone to 
have flux densities that are over-estimated due to stabilisation 
difficulties when the CVF scanned from the shortest towards 
longer wavelengths. 



B.2. PHOT-S spectro-photometry 

The measurements each consist of 4 or 8 DRs, with 127 NDRs 
per ramp, except for WOH G64 (16 DRs and 63 NDRs per 
ramp). All DRs were discarded, as well as the first 12 NDRs 
per ramp (6 in case of WOH G64) . Glitches were removed by 
applying the two-threshold median filtering technique, using 
distributions of 32 data points, three iterations, and thresholds 
of three standard deviations for flagging and re-accepting. The 
ramps were subdivided into four sections. The first four signals 
per measurement were discarded. The signals were deglitched 
by running a bin with a width of 16 signals over the chop- 
per plateau, taking individual steps, and iterating twice: af- 
ter maximum/minimum clipping, signals were discarded when 
they were off by more than three standard deviations. The sta- 
bility recognition method was applied, using an 8-signals wide 
bin and intervals of 4 signals, keeping only data within a con- 
fidence level of 0.95. The time evolution of the signal for each 
pixel was checked by eye, and obviously bad data that had 
passed the earlier rejection criteria were removed by hand at 
this stage. The expected dark current for the orbital position 
of the spacecraft was subtracted from the data. 

The in-band power was calibrated using the default respon- 
sivities expected for the orbital position, and finally the mean 
flux density was extracted, taking into account the amount of 
flux density lost outside of the aperture. Note that the wave- 
length dependence of the responsivity shows a bump around 3 
/jm. When applied to the spectrum of a star, this could intro- 
duce an artificial spectral feature which might be interpreted 
erroneously as absorption due, for instance, to H2O ice. 

The PHOT-S observations did not include a separate mea- 
surement of the background. We used the average of the 
background spectra measured in the CAM-CVF images (Fig. 
Bl, errobars indicate the standard deviation in the set of 
12 spectra). This spectrum was used for wavelengths > 8 
/xm, and interpolated linearly below 8 /im to a value of 
zero at 5 /im (at shorter wavelengths the background is as- 
sumed to be negligible). We ignore the first 4 steps of the 
CVF filter wheel, that are usually not stabilised. The back- 
ground is predominantly emission from the zodiacal dust 
belt in our Solar System. Because of the annual modula- 
tion of the zodiacal light level in a particular direction of 
the sky, the weekly all-sky maps produced by the DIRBE in- 
s trument onboard the COsmic Background Ex plorer (COBE) 
(tittp://www.gsfc.nasa.gov/astro/cobe/#dirbe) were used to 




scale the background spectrum to the epochs of the PHOT- 
S spectra (see Abraham et al. 1998). The COBE/DIRBE 12 
/im surface brightnesses for our CVF and PHOT-S observa- 



12 
X(fim) 

Fig. B.l. The CAM-CVF background emission spectrum. 
Open symbols represent spectro-photometric points that 
are prone to have flux densities that are over-estimated 
due to stabilisation difficulties. The spectral shape is best 
represented by the solid symbols (squares for the short- 
, disks for the long- wavelength region). The background 
emission actually arises from the zodiacal dust belt in our 
Solar system. 

tions have standard deviations of 5.5 and 5.9%, respectively. 
The CVF background spectra yield a standard deviation at 12 
/im of 9.1%. Although the CVF and COBE/DIRBE do corre- 
late, the spread in the CVF background exceeds the variation 
in the zodiacal light, suggesting that variations in the back- 
ground emission from the LMC is discernible in our data. The 
scaled background spectrum is subtracted from the PHOT-S 
spectra, that were obtained through a 24" x 24" aperture. 

IRAS05112— 6755 was observed twice. The spectra look 
identical, except that the second is a factor 1.5 brighter than 
the first, independent of wavelength, whilst the near-IR flux 
densities only differed by ~10%. Pointing thresholds were 10" 
and 2", respectively, and the star may have been close to 
the edge of the aperture during the first observation. Flux 
density levels in the PHOT-S spectrum of IRAS05298-6957 
may have been under-estimated due to off-centering of the 
star in the PHOT-S aperture. The PHOT-S spectrum of 
IRAS04496-6958 is fainter than the CAM-CVF spectrum, 
whereas for IRAS05558— 7000 the reverse is true. This may 
be a result of variability, as seen from the near-IR photometry 
at the epochs of these ISO observations. 
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